Hyperfine Field and Hyperfine Anomalies of Copper Impurities in Iron 
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A new value for the hyperfine magnetic field of copper impurities in iron is obtained by combining 
resonance frequencies from /3-NMR/ON experiments on ^^Cu, ''^Cu and ^^Cu with magnetic moment 
values from coUinear laser spectroscopy measurements on these isotopes. The resulting value, i.e. 
-Bhf(CuFe) = —21.794(10) T, is in agreement with the value adopted until now but is an order of 
magnitude more precise. It is consistent with predictions from ab initio calculations. Comparing 
the hyperfine field values obtained for the individual isotopes, the hyperfine anomalies in Fe were 
determined to be '^^A^^ = 0.15(9) % and "A"^^ = 0.07(11) %. 
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I. INTRODUCTION 

Precise values of magnetic hyperfine fields [1] allow the 
determination of nuclear magnetic moments by experi- 
mental methods such as integral perturbed angular cor- 
relation (IPAC) [2] and time-differential perturbed angu- 
lar distribution (TDPAD) [3] , or low-temperature nuclear 
orientation (LTNO) and nuclear magnetic resonance on 
oriented nuclei (NMR/ON) [4]. Precise results, in addi- 
tion, allow a detailed comparison with theory. 

The magnetic hyperfine fields of substitutional impu- 
rities in bcc Fe are at present well understood for most 
of the elements in the periodic table [5-11], with size- 
able differences between theory and experiment remain- 
ing mainly for the heavier 5d impurities [12, 13], the 
alkaline elements [14-17], and the actinides [11]. Still, 
for a few elements no precise experimental results are 
available yet. A special case is copper, for which the 
currently accepted value of the hyperfine field in iron is 
Shf(Cufe) = -21.8(1) T [18]. This value resulted from 
a spin-echo NMR (nuclear magnetic resonance) measure- 
ment with the sample containing the stable isotopes ^^Cu 
and ^^Cu cooled cooled to a temperature of 4.2 K. The er- 
ror of 0.1 T was not given in the original publication but 
was provided later as a private communication by one of 
the original authors [19]. The sign was obtained from the 
field shift in NMR measurements [20] and was confirmed 
by theoretical calculations (Ref. [21] and Sec. VI). 

In the past few years the magnetic hyperfine interac- 
tion frequencies z/ros ~ M^tot, with fj, the nuclear mag- 
netic moment and Btot the total magnetic field the nu- 
clei experience, have been determined for the Cu isotopes 
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^^Cu [22], 67Cu [19], 69Cu [23], and ^^Cu [24] with the 
^-NMR/ON method, i.e. NMR/ON with /3-particle de- 
tection. In these measurements, performed on samples 
that were cooled to millikelvin temperatures, the above 
mentioned value for the hyperfine field of Cu impurities 
in Fe host, viz. i3hf(Cui^e) = —21.8(1) T, was used in 
order to extract the nuclear magnetic moments for the 
isotopes studied. Recently, the magnetic moments of the 
Cu isotopes with A — 61 to 75, i.e. including the ones 
mentioned above, have been determined by coUinear laser 
spectroscopy measurements at ISOLDE-CERN [25-27]. 
As results for the /3-NMR/ON resonance frequencies and 
for the magnetic moments from the laser spectroscopy 
measurements for these isotopes all have similarly high 
precisions, ranging from 2 to 6 x 10~^, combining these 
results allows determining a new and precise value for the 
hyperfine magnetic field of Cu in Fe at K. 

Note that the resonance frequencies and magnetic mo- 
ment values that will be used here were all obtained with 
the same experimental methods and setups. This reduces 
the risk for possible systematic errors related to calibra- 
tion issues. 



II. NUCLEAR MAGNETIC RESONANCE ON 
ORIENTED NUCLEI 

The NMR/ON method has been appfied widely for the 
precise determination of the magnetic hyperfine splitting 
of radioactive nuclei in the ferromagnetic host lattices 
Fe, Co and Ni. In most cases, the primary goal was to 
deduce the nuclear magnetic moments of the impurity 
isotopes [28]. The same technique is also used to ground 
state spins [29], hyperfine fields [1, 30], nuclear relaxation 
times [31, 32], and quadrupole splittings [33, 34], as well 
as to obtain information on the lattice location and im- 
plantation behavior of implanted impurities [13, 35-37]. 



NMR/ON experiments require the nuclei to be oriented, 
which is done using the LTNO method [4] and requires 
cooling down the radioactive samples to temperatures in 
the millikelvin region and subjecting them to high mag- 
netic fields, either hyperfine magnetic fields [4] or exter- 
nally applied fields [38]. 

In NMR/ON, the resonant depolarization of the ra- 
dioactive probe nuclei is detected as a function of the 
frequency of the applied radio-frequency field via the re- 
sulting destruction of anisotropy in the anisotropic emis- 
sion of the decay radiation. The resonance frequency is 
related to the hyperfine field, i?hf , through the relation 



J^rcs[MHz] 



7.6226 -AiM^totiTesla] 



m 



with I the nuclear spin of the isotope studied and 



Btot ^ Bm + B,pp{l + K) - B. 



^dcm • 



(1) 



(2) 



Here the hyperfine field B^ includes the Lorentz field of 
0.742 T for bcc iron at K, i?app is the externally applied 
magnetic field, i?dem is the demagnetization field and, K 
is the Knight shift. The constant factor is the ratio of the 
fundamental constants ^n /h and does not contribute to 
the error budgets here. 



III. COLLINEAR LASER SPECTROSCOPY 

CoUinear laser spectroscopy experiments determine the 
hyperfine structure of atomic ground- and excited states, 
yielding precise values for the hyperfine parameters A 
and B, which in turn provide accurate values for mag- 
netic and quadrupole moments of the isotopes stud- 
ied [41]. 

Recently, collinear laser spectroscopy experiments were 
performed on a series of Cu isotopes at ISOLDE with 
the COLLAPS setup, including the four isotopes for 
which resonance frequencies are available from recent j3- 
NMR/ON measurements. For several of these isotopes 
laser spectroscopy could only be performed after the 
installation of the ISCOOL cooler and buncher radio- 
frequency quadrupole Paul trap [42, 43]. This allowed 
the collinear laser setup to be operated in bunched mode, 
reducing the background photon counts from scattered 
laser light by more than three orders of magnitude, and 
permitting to perform measurements on isotopes that 
were previously not accessible due to their low yields. 
The Cu"*" bunches from ISCOOL were sent through a 
sodium vapor cell which neutralized the ions through 
charge- exchange collisions. A voltage was applied to the 
vapor cell for tuning the velocity of the ions and bringing 
them onto resonance with the laser beam that was over- 
lapped with the Cu beam in the co-propagating direc- 
tion. Resonances were located by measuring the photon 
yield as a function of the voltage with two photomulti- 
plier tubes, the voltage of which was gated so that pho- 
tons were only recorded when an atom bunch was within 



the light collection region [25] . An example of a collinear 
resonance fluorescence spectrum obtained in these mea- 
surements is shown in Fig. 1 of Rcf. [25] . 



IV. HYPERFINE MAGNETIC FIELD FOR 
COPPER IN IRON 

In the following the magnetic moments from the 
collinear laser spectroscopy measurements on the iso- 
topes 59Cu [27], 6^Cu [26], s^Cu [26] and ^^Cu [25] are 
combined with the /3-NMR/ON measurements that were 
performed on these isotopes [19, 22-24]. The experimen- 
tal magnetic moment values and /3-NMR/ON resonance 
frequencies for ^^Cu, ^®Cu, and ^^Cu are listed in Table I 
and will be discussed in Sec. IV B. The data for ^^Cu are 
discussed separately in Sec. IV C. 

As the extraction of the hyperfine magnetic field using 
Eqs. 1 and 2 requires that a possible Knight shift and the 
demagnetization field are taken into account as well, the 
/3-NMR/ON measurements and these two factors will be 
addressed in some detail first. 



A. /3-NMR/ON measurements 

All four Cu isotopes on which /3-NMR/ON mea- 
surements have been performed were produced at the 
ISOLDE isotope separator facility The NMR/ON 
measurements were performed either on-line using the 
NICOLE LTNO setup [44] at ISOLDE (for ^^Cu [22], 
69 Cu [23] and ^^Cu) [24], or off-line with the nuclear- 
orientation facility at Oxford University (for ^^Cu [19]). 
As an example, the NMR/ON result for ^^CuFe is shown 
in Fig. 1. 



1. Demagnetization field 

The calculation of the demagnetization field is not 
straightforward, but for the simple shapes and thin foils 
(magnetized in the plane of the foil) that were used in 
the /3-NMR/ON measurements discussed here, analyti- 
cal expressions can be obtained [40] . At the current level 
of precision, the small demagnetization fields in the foils 
used, which are of the order of 0.01 T, cannot be ne- 
glected. 

The foil used by the Lcuven group for the measure- 
ment with ^^Cu had a size of 9 mm by 14 mm and 
an initial thickness of 125 /xm. It was polished with 
diamond-base paste with grain sizes of 3 /im and 1 jjlth. 
It is estimated that this reduced the thickness of this 
foil to (95 ± 20) /im. For this a demagnetization field 
-Bdcm = 0.020(4) T is calculated. The foils used by the 
Oxford group for the measurements with ^ Cu, ^^Cu and 
^^Cu had an initial thickness of 250 /xm [45]. They were 
first polished with sand paper with CAMI (Coated Abra- 
sive Manufacturers Institute) grit designation 600 (i.e. 
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FIG. 1. On-line NMR/ON curve for ^^Cu isotope. Plotted 
is the ratio of the count rates for two detectors at 15° and 
165° with respect to the quantization axis as a function of 
frequency. At the bottom, the count rate at 0° is shown for 
the 136 keV 7 ray of the ^^CoFe thermometer (corresponding 
to a sample temperature of about 10 mK) for the same fre- 
quency region, showing no effect at the position of the *^Cu 
resonance. The slope in the count rates is caused by a small 
increase of the sample temperature with the resonance fre- 
quency, (arrows indicate the vertical scale that applies; from 
ref. [22]) 



with an average particle diameter of 16 /im) and there- 
after also with diamond-base paste. Due to the larger 
particle diameter used in the first step, this procedure 
removed more material from the foils leading to an esti- 
mated thickness of (190 ± 40) \iva.. For typical dimen- 
sions of 10 mm by 15 mm for such foils, a demagnetiza- 
tion field Bdom = 0.036(8) T is then calculated for the 
measurements with ^^Cu, ^^Cu and ^^Cu. 



2. Knight shift 



The Knight shift for copper in iron has never been 
determined at low temperatures. For other elements in 
ferromagnetic host materials Knight shift values ranging 
from zero to about 5% have been reported [39]. Thus, for 
the low external fields, i?app, of 0.1 T and 0.2 T that were 
used in the /3-NMR/ON measurements discussed here, 
the Knight shift corrections could amount up to about 
0.005 T and 0.010 T, respectively. These values were then 
used as the la error bars on the values for the externally 
applied magnetic fields. Note that the precision to which 
these external fields could be set is an order of magnitude 
better and can therefore be neglected at the current level 
of precision. 



B. Results for '^^Cu, ®^Cu and "Cu 

The magnetic moment values and the resonance fre- 
quencies for the isotopes ^^Cu, ^^Cu and ^^Cu are listed 
in Table I (the case of ^^Cu will be discussed separately 
in the next section). Also listed there are the values for 
the total magnetic field, i?tot, obtained for each isotope 
using Eq. 1, as well as the hyperfine field, -Bhf, that is 
then obtained from Eq. 2 using the values for the exter- 
nally applied field, -Bapp, and the demagnetization field, 

-^dcin- 

As can be seen, the hyperfine field values obtained for 
these three isotopes are in very good agreement with each 
other (see also Fig. 2). Combining all three results yields 
a weighted average value of 



5hf(CuFe) 



-21.794(10) T 



(3) 



The differences between the hyperfine field values ob- 
tained for the three isotopes listed in Table I could be 
due to small differences in the distribution of the nuclear 
magnetization over the nuclear volume for the three iso- 
topes, i.e. hyperfine anomalies [46]. The usual definition 
for the hyperfine anomaly, e^, for a single nuclear state i 



B, 



off 



Bo(l 



(4) 



where Bes is the hyperfine field averaged over the distri- 
bution of nuclear magnetization of the state i and Bq is 
the hyperfine field at r = 0. The difference of the hy- 
perfine anomalies of two nuclear states in the same host 
metal is then given by 



^A^ 



ei - £2 






1. 



(5) 



When considering the possible presence of hyperfine 
anomalies the values for Bm listed in Table I have to be 
interpreted as values for i?cff. It then turns out that the 
differences of the hyperfine anomalies of these three iso- 
topes are less than 3 x 10"^ (see Table I; 90% C.L.). Note 
that previously the hyperfine anomaly between ^'^Cu and 
^^Cu was estimated to be less than 5 x 10~^ [47], while 
in the recent coUinear laser spectroscopy measurements 
on Cu isotopes [26], no indication for hyperfine anoma- 
lies was observed either for the series of isotopes ranging 
from A = 61 to 75. 

For ^^Cu, two other values for the mag- 
netic moment were recently quoted as well, i.e. 
/i(^^Cu) = 1.83(4) fiN, from in-source laser spectroscopy 
[48], and /i(^^Cu) = 1.910(4) ^n, from an in-gas-cell 
laser spectroscopy experiment [49, 50]. As the first value 
is much less precise than the one obtained from collinear 
laser spectroscopy it is not further used here. The second 
value differs well outside error bars from the collinear 
laser spectroscopy result (see Table I). Combining it with 
the /3-NMR/ON resonance frequency of 208.79(4) MHz 
for ^^Cu reported in Ref. [22], a hyperfine field value 



of Bhf(CuFe) = -21.59(5) T is obtained. This differs 
by about 4 standard deviations from the values hsted 
in Table I that resulted from combining resonance 
frequencies and magnetic moment values that were 
each obtained with the same experimental setups. This 
difference might be the due to unforeseen systematic 
effects in this measurement [49] which was performed 
with a different experimental setup, as is currently being 
investigated [51]. 



C. The case of "Cu 

The situation for ^^Cu turns out to be more com- 
plicated. In Ref. [23], the magnetic moment of ^^Cu 
obtained with the Oxford LTNO setup is given as 
/^(67Cu) = +2.54(2) /xjv; in Ref. [52] this result is quoted 
as /x = +2.536(3) fiN (this smaller error bar most prob- 
ably only reflects the statistical precision). The authors 
do not quote a resonance frequency value, but analysis 
of the resonance curve shown in Ref. [23] yields a central 
value of t'ros = 278.38(6) MHz. In combination with the 
magnetic moment value of /i = +2.536 fiN, this yields a 
total magnetic field of —21.60 T. Comparing this with the 
hyperfine field value of —21.8(1) T for copper in iron that 
was used in Ref. [23] and neglecting the demagnetization 
field (as the authors of Ref. [23] also did), yields a value of 
+0.20 T for the externally applied magnetic field instead 
of the value of +0.10 T mentioned in Refs. [23, 52]. The 
latter value is therefore most probably a typographical 
error. 

When combining the above resonance frequency of 
278.38(6) MHz with the magnetic moment value 
/i(^^Cu) — +2.5142(6) fiN from the laser spectroscopy 
experiments [26], a hyperfine field value of i?hf (CuFe) = 
-21.952(15) T is obtained for Boxt = 0.20(1) T and with 
-Sdom = 0.036(8) T for the foils of the Oxford team (see 
Sec. IVAl). If the demagnetization field is neglected 
Bhf(CuFe) = -21.988(15) T is obtained. Both values 
are at variance with the ones obtained for the other Cu 
isotopes (see Table I and Fig. 2). 

Concluding, there seems to be a problem with the /3- 
NMR/ON result for ^^Cu that was obtained using an- 
other experimental setup than the one used for the iso- 
topes ^^Cu, ^^Cu and ^^Cu. The origin of this may e.g. 
be an undetected error in the frequency calibration. We 
therefore did not include the /?-NMR/ON result for •'^Cu 
in the hyperfine field analysis presented here. 



V. PREVIOUS RESULTS 

Table II lists other values for the hyperfine field of 
Cu in Fe that are available in the literature, most of 
which were obtained in classical NMR experiments at 
room temperature and with - unfortunately - no error 
being quoted. As can be seen, most values are in reason- 
able agreement with the new value of —21.794(10) T pre- 
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FIG. 2. Hyperfine field values (in Tesia) for Cu in Fe ob- 
tained by combining resonance frequencies from /3-NMR/ON 
measurements on the isotopes ^^Cu, ®^Cu, ''^Cu and ^^Cu 
with the magnetic moments for these isotopes from collinear 
laser spectroscopy. The grey band indicates the weighted av- 
erage of the values from Cu, 



sented here. The value from the 7-7 perturbed angular 
correlation measurement, which is significantly deviating 
from all other results, is either wrong or might be related 
to a different lattice site for Cu impurities in Fe. 



VI. COMPARISON WITH THEORETICAL 
VALUES 

Hyperfine fields in solids can be calculated from first 
principles. This allows understanding trends in those hy- 
perfine fields as a function of the impurity element in e.g. 
an iron matrix [5-11], and it allows disentangling the hy- 
perfine field into contributions with a different physical 
origin [10, 11, 61, 62]. We have calculated the hyperfine 
field of Cu in Fe within the framework of Density Func- 
tional Theory [63-65] , using the Perdew-Burke-Ernzerhof 
(FEE) exchange-correlation functional [66]. For solving 
the scalar-relativistic Kohn-Sham equations we employed 
the Augmented Plane Waves + local orbitals (APW+lo) 
method [65, 67, 68] as implemented in the wiEN2k pack- 
age [69] for periodic solids. In this method the wave 
functions are expanded into spherical harmonics inside 
nonoverlapping atomic spheres of radius -Rmt and in 
plane waves in the remaining space of the unit cell, i.e. 
the interstitial region. We took i?MT — -^mt = 2.30 a.u. . 
The plane wave expansion of the wave function in the in- 
terstitial region was truncated at a large value of A'max = 
8.0/i?™^ — 3.48 a.u.~^, which leads to very accurate val- 
ues for the calculated hyperfine fields. A dense mesh of 
k-points, corresponding to a 20x20x20 mesh for a con- 
ventional bcc unit cell for Fe, was taken. Spin-orbit cou- 
pling was taken into account by a second variational step 
scheme [70] using a cutoff energy Eso= 5.0 Ry. The 



TABLE I. Data used for the determination of the hyperfine field of Cu impurities in Fe and values deduced from these. The 
error bar on the weighted average was increased by a factor ^y} jv = \/1.22 to account for the fact that the reduced y^ is 
larger than unity. The last column lists the hyperfine anomaly differences between ^^Cu and the two other isotopes. 



Isotope 


(^iv) 


Rcf. 


(MHz) 


Rcf. 


(T) 


R b 

^app 

(T) 


Sdom 

(T) 


(T) 


(%) 


s^Cu 
e^Cu 
"Cu 


+1.8910(9) 
+2.8383(10) 
+2.2747(8) 


[27] 

[26] 
[25] 


208.79(4) 

311.7(1) 

250.00(14) 


[22] 

[23] 
[24] 


-21.727(11) 
-21.611(10) 
-21.627(14) 


0.100(5) 
0.20(1) 
0.20(1) 


0.020(4)^= 
0.036(8)" 
0.036(8)'= 


-21.807(13) 
-21.775(16) 
-21.791(19) 


0.15(9) 
0.07(11) 


Weighted average 
















-21.794(10) 





^Sign from [57]. '^From the Refs. listed in column 3. '^Calculated with the formulas derived in Ref. [40] (see also Sect. IV A 1) 



TABLE II. Overview of hyperfine magnetic fields for Cu in 
Fe host reported in the literature. Note that in most cases no 
error bar was quoted. 



Bhf (T) 


Rcf. 


method 


21.0 


[53] 


NMR 


21.77 


[54] 


NMR 


-21.27 


[55] 


NMR 


21.79 


[56] 


NMR 


-21.3 


[57] 


spin-echo 


21.8(1) 


[18, 19] 


NMR 


21.3 


[58] 


NMR 


-21.84 


[59] 


NMR 


16.95(51) 


[60] 


7-7 PAC" 



^perturbed angular correlation. 

substitutional Cu impurity was modeled by a 128- atom, 
supercell, and all atoms in the supercell were allowed to 
adjust their positions due to the presence of the impu- 
rity. The lattice constant for the Fe-matrix was taken 
to be the equilibrium lattice constant for the PBE func- 
tional (2.8404 A). These settings allow for an excellent 
numerical convergence of the hyperfine field. 

The results of the calculations are summarized in Ta- 
ble III. As can be seen, the distance between the Cu 



TABLE III. Calculated values for the Fermi contact field, or- 
bital hyperfine field, dipolar hyperfine field and atomic dipolar 
hyperfine field for a Cu impurity in Fe. Spin and orbital con- 
tributions to the atomic magnetic moment for Cu (/is) and 
distance between Cu and its first Fe-neighbor (A). 



Bf 


-25.32 


/^spin 


0.115 


Borb 


0.60 


/iorb 


0.007 


Bdip 


-0.01 


rfCu-Fe 


2.478 


BhL 


0.00 







electron) magnetic moment at the Cu atom. This value 
can be further split into a core contribution due to Is and 
2s electrons (—8.54 T), a semi-core contribution by the 
3s-clectrons (+5.80 T) and a valence contribution by 4s 
electrons (—22.59 T). The result for the Fermi contact hy- 
perfine field can be compared with the —18.2 T that was 
obtained 25 years ago by the Korringa-Kohn-Rostoker 
Green's function method [6]. The orbital hyperfine field 
of Cu (0.60 T) is an order of magnitude smaller than the 
corresponding quantity in pure Fe, consistent with the 
very small atomic orbital magnetic moment of Cu. All 
contributions to the hyperfine field, together with the 
Lorenz field (0.74 T), sum to a total hyperfine field of 
—23.99 T. Although this is about 2 T larger than the ex- 
perimental value determined in this work, this deviation 
is state-of-the-art, and is due to inherent limitations of 
the chosen exchange-correlation functional. 



VII. CONCLUSION 

Combining resonance frequencies for the isotopes ^^Cu, 
^^Cu and ^^Cu obtained from /3-NMR/ON measure- 
ments with the NICOLE LTNO setup at ISOLDE, with 
magnetic moment values obtained for these isotopes in 
coUinear laser spectroscopy measurements at ISOLDE, 
the hyperfine field of Cu impurities in iron is found to be 
Shf (CuFe) = -21.794(10) T. This value is in agreement 
with but almost an order of magnitude more precise than 
the previously adopted value of —21.8(1) T and in good 
agreement with predictions from ah initio calculations. 
Interpreting the differences between the hyperfine field 
values obtained for the individual isotopes to be due to 
hyperfine anomalies, the hyperfine anomalies in Fe for the 
isotopes considered here were found to be smaller than 
3 X 10-3 (90% C.L.). 



impurity and its first eight Fe neighbours is expanded 
only slightly (0.75 %) compared to the Fe-Fe distance 
of 2.460 A in pure Fe. The dominant contribution to 
the hyperfine field is the Fermi contact term, caused by 
s-electron spin polarization, due to the small atomic (d- 



VIII. ACKNOWLEDGEMENT 

This work was supported by the Fund for Scientific 
Research Flanders (FWO), project GOA/2004/03 of the 
K. U. Leuven, the Interuniversity Attraction Poles Pro- 



gramme, Belgian State Belgian Seience Poliey (BriX net- 



work P6/23), and the grant LA08015 of the Ministry of 
Edueation of the Czeeh Republic. 



[1] G. Rao, Hyperfine Interact. 24-26, 1119 (1985). 
[2] E. Bodenstedt, in The Electromagnetic Interaction in 
Nuclear Spectroscopy, edited by W.D. Hamilton (North- 
Holland, Amsterdam, 1975), p. 735. 
[3] G. Goldring and M. Hass, in Treaties on Heavy Ion Sci- 
ence Vol. 3, edited by D. Allan Bromley (Plenum, New 
York, 1985) 539. 
[4] N. J. Stone in Low-Temperature Nuclear Orientation, ed- 
itors H. Postma and N. J. Stone (North-Holland, Ams- 
terdam, 1986). 
[5] H. Akai, M. Akai, S. Bliigel, R. Zeller, and P.H. Ded- 

erichs, J. Phys. Soc. Jpn. 45, 291 (1984). 
[6] M. Akai, H. Akai, and J. Kanamori, J. Phys. Soc. Jpn. 

54, 4246 (1985). 
[7] M. Akai, H. Akai, and J. Kanamori, J. Phys. Soc. Jpn. 

54, 4257 (1984). 
[8] T. Korhonen, A. Settels, N. Papanikolaou, R. Zeller, and 

P.H. Dederichs, Phys. Rev. B 62, 452 (2000). 
[9] S. Cottenier and H. Haas, Phys. Rev. B 62, 461 (2000). 

[10] D. Torumba, V. Vanhoof, M. Rots, and S. Cottenier, 
Phys. Rev. B 74, 014409 (2006). 

[11] D. Torumba, P. Novak, and S. Cottenier, Phys. Rev. B 
77, 155101 (2008). 

[12] H. Ebert, R. Zeller, B. Drittler, and P.H. Dederichs, J. 
Appl. Phys. 67, 4576 (1990). 

[13] N. Severijns et al, Phys. Rev. C 79, 064322 (2009). 

[14] J. Wouters, N. Severijns, D. Vandeplassche, E. vanWalle, 
and L. Vanneste, Phys. Lett. A 124, 377 (1987). 

[15] W. Vanderpoorten, P. De Moor, P. Schuuramns, R. 
Siebelink, L. Vanneste, J. Wouters, N. Severijns, J. Van- 
haverbeke, R. Eder, and H. Haas, Hyperfine Interact. 75, 
331 (1992). 

[16] C. J. Ashworth, P. Back, S. Ohya, N. J. Stone and J. P. 
White, Hyperfine Interact. 59, 461 (1990). 

[17] B. Will et al, Phys. Rev. B 57, 11527 (1998). 

[18] L. Khoi, P. Veillet, and I. CampbeU, J. Phys. F 5, 2184 
(1975). 

[19] J. Rikovskaand N. J. Stone, Hyperfine Interact. 129, 131 
(2000). 

[20] M. Kontani and J. Itoh, J. Phys. Soc. Japan 22, 345 
(1967). 

[21] H. Akai, M. Akai, and J. Kanamori, J. Phys. Soc. Jpn 
54, 4257 (1985). 

[22] V. V. Golovko et al, Phys. Rev. C 70, 014312 (2004). 

[23] J. Rikovska et al, Phys. Rev. Lett. 85, 1392 (2000). 

[24] N. J. Stone et al., Phys. Rev. C 77, 014315 (2008). 

[25] K. T. Flanagan et al, Phys. Rev. Lett. 103, 142501 
(2009). 

[26] P. Vingerhoets et al, Phys. Rev. C 82, 064311 (2010). 

[27] P. Vingerhoets, PhD thesis, Kath. University Leuven 
(2011), and submitted to Phys. Lett. A. 

[28] N.J. Stone, Atomic Data and Nuclear Data Tables 90, 
75 (2005). 

[29] D. Vandeplassche, E. van Walle, J. Wouters, N. Severijns 
and L. Vanneste, Phys. Rev. Lett. 57, 2641 (1986). 

[30] K. S. Krane, Hyperfine Interact. 15/16, 1069 (1983). 



[31] F. Bacon, J. A. Barclay, W. D. Brewer, D. A. Shirley, 

and J. E. Templeton, Phys. Rev. B 5, 2397 (1972). 
[32] P. Herzog in Ref. [4], p. 953. 
[33] P. D. Johnston, R. A. Fox, and N. J. Stone, J. Phys. C, 

5, 2077 (1972). 
[34] E. Hagn, Hyperfine Interact. 22, 19 (1985). 
[35] H. Pattyn et al, Hyp. Int. 2, 362 (1976). 
[36] U. Dammrich and P. Herzog, Hyperfine Interact. 43, 169 

(1988). 
[37] E. van Walle, D. Vandeplassche, J. Wouters, N. Severijns, 

and L. Vanneste, Phys. Rev. B 34, 2014 (1986). 
[38] W. D. Brewer in Ref. [4], Chap. 9. 
[39] N. Yazidjoglou, W. D. Hutchison, and D. H. Chaplin, J. 

Phys. Condens. Matter 5, 129 (1993). 
[40] S. Chikazumi, Physics of magnetism, Vol. 1 (J. Wiley and 

Sons, New York, 1964). 
[41] A. Mueller et al., Nucl. Phys. A 403, 234 (1983). 
[42] H. Franberg et al, Nucl. lustrum. Methods Phys. Res. B 

204, 4502 (2008). 
[43] E. Mane et al., Eur. Phys. J. A 42, 503 (2009). 
[44] K. Schlosser, I. Berkes, E. hagn, P. Herzog, T. Niinikoski, 

H. Postma, C. Richard- Serre, J. Rikovska, N. J. Stone, 

L. Vanneste, E. Zech and the NICOLE and ISOLDE Col- 
laborations, Hyperfine Interact. 43, 141 (1988). 
[45] K. Van Esbroeck, PhD thesis, Cath. University of Leuven 

(2005), unpublished. 
[46] S. Biittgenbach, Hyperfine Interact. 20, 1 (1984). 
[47] P. R. Locher, Phys. Rev. B 10, 801 (1974). 
[48] N. J. Stone, U. Koster, J. R. Stone, D. V. Fedorov, 

V. N.Fedoseyev, K. T. Flanagan, M. Hass, and S. Lak- 

shmi, Phys. Rev. C 77, 067302 (2008). 
[49] T. E. Cocolios et al, Phys. Rev. Lett. 103, 102501 

(2009). 
[50] T. E. Cocolios et al., Phys. Rev. C 81, 014314 (2010). 
[51] T. E. Cocolios, private communication (2011). 
[52] J. Rikovska et al, CERN document CERN/ISC 99-6, 

ISC/P-79 Addendum 1. 
[53] Y. Koi, T. Kushida, T. Hihara, and A. Tsujimura, J. 

Phys. Soc. Jpn. 17, 96 (1962). 
[54] T. Kushida, A. H. Silver, Y. Koi, and A. Tsujimura, J. 

Appl. Phys. 33, 1079 (1962). 
[55] D. A. Shirley and G. A. Westenbarger, Phys. Rev. 138, 

A170 (1965). 
[56] D. Edmonds and G. Wilson, Phys. Lett. 23, 431 (1966). 
[57] M. Kontani and J. Itoh, J. Phys. Soc. Jpn. 22, 345 (1967). 
[58] P. Ricdi and G. Webber, J. Phys F 11, 1669 (1981). 
[59] Y. Kasamatsu, T. Hihara, K. Kojima, and T. Kami- 

gaichi, J. Magn. Magn. Mater. 54-57, 1107 (1986). 
[60] E. Lohmann, K. Freitag, T. Schaefer, and R. Vianden, 

Hyperfine Interact. 77, 103 (1993). 
[61] P. Novak, J. Kunes, W.E. Pickett, Wei Ku, and F. R. 

Wagner, Phys. Rev. B 67, 140403(R) (2003). 
[62] E. L. Peltzer y Blanca, J. Desimoni, N. E. Christensen, 

H. Emmerich, and S. Cottenier, Phys. Status Solidi B 

246, 909 (2009). 
[63] P. Hohonborg and W. Kohn, Phys. Rev. 136, B864 

(1964). 



[64] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

[65] S. Cottenier, Density Functional Theory and the Fam- 
ily of (L)APW-methods: a step-by-step introduction, In- 
stituut voor Kern- en Stralingsfysica, KULeuven, Bel- 
gium, ISBN 90-807215-1-4 (freely available from http: 
//www . wien2k . at/reg_user /textbooks) . 

[66] J. P. Perdew, K. Burke, and M. Ernzcrhof, Physical Re- 
view Letters 77, 3865 (1996). 

[67] E. Sjostedt , L. Nordstrom, and D. J. Singh, Solid State 
Communications 114, 15 (2000). 



[68] G. K. H. Madsen, P. Blaha, K. Schwarz, E. Sjostedt, and 
L. Nordstrom, Physical Review B 64, 195134 (2001). 

[69] P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. 
Luitz, WIEN2k, An Augmented Plane Wave-|-Local Or- 
bitals Program for Calculating Crystal Properties (Karl- 
heinz Schwarz, Technische Universitt Wien, Austria, 
1999). 

[70] D. D. Koelling and B. N. Harmon, J. Phys. C 10, 3107 
(1977). 



